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Abstract 

We consider the problem of cutting a set of edges on a polyhedral manifold surface, 
possibly with boundary, to obtain a single topological disk, minimizing either the total 
number of cut edges or their total length. We show that this problem is NP-hard, 
even for manifolds without boundary and for punctured spheres. We also describe an 
algorithm with running time n°^ 9+k \ where n is the combinatorial complexity, g is the 
genus, and k is the number of boundary components of the input surface. Finally, we 
describe a greedy algorithm that outputs a 0(log 2 (^-approximation of the minimum 
cut graph in 0(g 2 n log n) time. 
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1 Introduction 

Several applications of three-dimensional surfaces require information about the underlying 
topological structure in addition to the geometry. In some cases, we wish to simplify the 
surface topology to facilitate algorithms that can be performed only if the surface is a 
topological disk. 

Applications when this is important include surface parameterization [20, 46] and texture 
mapping [3, 44]. In the texture mapping problem, we wish to find a continuous and invertible 
mapping from the texture, usually a two-dimensional rectangular image, to the surface. 
Unfortunately, if the surface is not a topological disk, no such map exists. In such a case, 
the only feasible solution is to cut the surface so that it becomes a topological disk. (Haker 
et al. [26] present an algorithm for directly texture mapping models with the topology of a 
sphere, where the texture is also embedded on a sphere.) Of course, when cutting the surface, 
one would like to find the best possible cut under various considerations. For example, one 
might want to cut the surface so that the resulting surface can be textured mapped with 
minimum distortion [20, 25, 46]. To our knowledge, all previous approaches for this cutting 
problem either rely on heuristics with no quality guarantees [25, 45, 48] or require the user 
to perform this cutting beforehand [20, 44]. 

Lazarus et al. [37] presented and implemented two algorithms for computing a canonical 
polygonal schema of an orientable surface of complexity n and with genus g, in time O(gri), 
simplifying an earlier algorithm of Vegter and Yap [55]. Computing such a schema requires 
finding 2g cycles, all passing through a common basepoint in A4, such that cutting along 
those cycles breaks M. into a topological disk. Since these cycles must share a common point, 
it is easy to find examples where the overall size of those cycles is Vt(gn). Furthermore, those 
cycles share several edges and are visually unsatisfying. 

For most applications, computing a canonical schema is overkill. It is usually sufficient 
to find a collection of edges whose removal transforms the surface into a topological disk. 
We call such a set of edges a cut graph; see Figure 1 for an example. Cut graphs have several 
advantages. First, they are compact. Trivially, any cut graph contains at most n edges of 
the surface mesh, much less than any canonical schema in the worst case, although we expect 
it to be much smaller in practice. Second, it is quite easy to construct a cut graph for an 
arbitrary polyhedral surface in 0{n) time, using a breadth-first search of the dual graph [13], 
or simply taking a maximal set of edges whose complement is connected [37]. Finally, the 
cut graph has an extremely simple structure: a tree with 0(g) additional edges. As such, 
it should be easier to manipulate algorithmically than other representations. For example, 
Dey and Schipper [13] describe fast algorithms to determine whether a curve is contractible, 
or two curves are homotopic, using an arbitrary cut graph instead of a canonical schema. 

In this paper, we investigate the question of how find the "best" such cutting of a surface, 
restricting ourselves to cuts along the edges of the given mesh. Specifically, we want to find 
the smallest subset of edges of a polyhedral manifold surface Ai, possibly with boundary, 
such that cutting along those edges transforms M. into a topological disk. We also consider 
the weighted version of this problem, where each edge has an arbitrary non-negative weight 
and we want to minimize the total weight of the cut graph. The most natural weight of 
an edge is its Euclidean length, but we could also assign weights to take problem-specific 
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Figure 1. A cut graph for a two-holed torus and its induced (non-canonical) polygonal schema. 

considerations into account. For example, if we want to compute a texture mapping for a 
specific viewpoint, we could make visible edges more expensive, so that the minimum cut 
graph would minimize the number of visible edges used in the cuts. Our algorithms do not 
require the edge weights to satisfy the triangle inequality. 

We show that the minimum cut graph of any polyhedral manifold M. with genus g and 
k boundary components can be computed in n 0( - 9+kS> time. We also show that the problem is 
NP-hard in general, even if g or k is fixed. Finally, we present a simple and efficient greedy 
approximation algorithm for this problem. Our algorithm outputs a cut graph whose weight 
is a factor 0(log 2 g) larger than optimal, in 0(g 2 n log n) time. 1 If g = 0, the approximation 
factor is exactly 2. As a tool in our approximation algorithm, we also describe efficient 
algorithms to compute shortest and nearly-shortest nontrivial cycles in a manifold; we believe 
these algorithms are of independent interest. 

2 Background 

Before presenting our new results, we review several useful notions from topology and 
describe related results in more detail. We refer the interested reader to Hatcher [29], 
Munkres [43], or Stillwell [49] for further topological background and more formal defini- 
tions. For related computational results, see the recent surveys by Dey, Edelsbrunner, and 
Guha [11] and Vegter [54]. 

2.1 Topology 

A 2-manifold with boundary is a set M. such that every point x G M. lies in a neighborhood 
homeomorphic to either the plane IR 2 or a closed halfplane. The points with only halfplane 
neighborhoods constitute the boundary of Ai; the boundary consists of zero or more disjoint 
circles. This paper will consider only compact manifolds, where every infinite sequence of 
points has a convergent subsequence. 

1 To simplify notation, we define log a; = max{l, |~log 2 x~\}. 
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The genus of a 2-manifold Ai is the maximum number of disjoint non-separating cycles 
71, 72, . . . , 7 9 in Ai; that is, 7$ fl 7 .,• = for all % and j, and Ai \ (71 U • • • U 7 9 ) is connected. 
For example, a sphere and a disc have genus 0, a torus and a Mobius strip have genus 1, and 
a Klein bottle has genus 2. 

A manifold is orientable if it has two distinct sides, and non-orientable if it has only 
one side. Although many geometric applications use only orientable 2-manifolds (primarily 
because non-orientable manifolds without boundary cannot be embedded in IR 3 without self- 
intersections) our results will apply to non-orientable manifolds as well. Every (compact, 
connected) 2-manifold with boundary is characterized by its orientability, its genus g, and 
the number k of boundary components [21]. 

A polyhedral 2-manifold is constructed by gluing closed simple polygons edge-to-edge 
into a cell complex: the intersection of any two polygons is either empty, a vertex of both, 
or an edge of both. We refer to the component polygons has facets. (Since the facets are 
closed, every polyhedral manifold is compact.) For any polyhedral manifold Ai, the number 
of vertices and facets, minus the number of edges, is the Euler characteristic x of Ai . Euler's 
formula [18] implies that x is an invariant of the underlying manifold, independent of any 
particular polyhedral representation; x = 2 — 2g — k if the manifold is orientable, and 
X = 2 — g — k if the manifold is non-orientable. Euler's formula implies that if Ai has v 
vertices, then Ai has at most 3v — 6 + Qg edges and at most 2v — 4 + Ag — k facets, with 
equality for orientable manifolds where every facet and boundary circle is a triangle. We let 
n < 6v — 10 + lOg — k denote the total number of facets, edges, and vertices in Ai. 

The 1-skeleton Ati of a polyhedral manifold Ai is the graph consisting of its vertices and 
edges. We define a cut graph G of Ai as a subgraph of Ai\ such that Ai \G is homeomorphic 
to a disk. 2 The disk Ai \ G is known as a polygonal schema of Ai. Each edge of G appears 
twice on the boundary of polygonal schema Ai \ G, and we can obtain Ai by gluing together 
these corresponding boundary edges. Finding a cut graph of Ai with minimum total length 
is clearly equivalent to to finding a polygonal schema of Ai with minimum perimeter. 

Any 2-manifold has a so-called canonical polygonal schema, whose combinatorial struc- 
ture depends only on the genus g, the number of boundary components k, and whether the 
manifold is orientable. 3 The canonical schema of an orientable manifold is a (4g + 3/c)-gon 
with successive edges labeled 

%1, 2/1 5 %1, Vli ■ ■ ■ 1 %g, ygi %g, ygi Zl, Zl, ■ ■ ■ ; ^fc; Cfc; ^fcj 

for a non-orientable manifold, the canonical schema is a (2g + 3/c)-gon with edge labels 

•^lj • • • 5 ^gi "Egi ^1) ^1) • • • j Z>ki Cfc) %k- 

Every pair of corresponding edges x and x is oriented in opposite directions. Gluing to- 
gether corresponding pairs in the indicated directions recovers the original manifold, with 

2 Cut graphs are generalizations of the cut locus of a manifold M., which is essentially the geodesic medial 
axis of a single point. 

3 Actually, there are several different ways to define canonical schemata; the one described here is merely 
the most common. For example, the canonical schema for an oriented surface without boundary could also 
be labeled x\,X2, ■ • ■ , X2 g , ah, x 2 , . . . , X2 3 - 
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the unmatched edges forming the boundary circles. For a manifold M. without bound- 
ary, a reduced polygonal schema is one where all the vertices are glued into a single point 
in Ai; canonical schemata of manifolds without boundary are reduced. We emphasize that 
the polygonal schemata constructed by our algorithms are neither necessarily canonical nor 
necessarily reduced. 

2.2 Previous and Related Results 

Dey and Schipper [13] describe an algorithm to construct a reduced, but not necessarily 
canonical, polygonal schema for any triangulated orientable manifold without boundary in 
0(n) time. Essentially, their algorithm constructs an arbitrary cut graph G by depth-first 
search, and and then shrinks a spanning tree of G to a single point. (See also Dey and 
Guha [12].) 

Vegter and Yap [55] developed an algorithm to construct a canonical schema in optimal 
O(gn) time and space. Two simpler algorithms with the same running time were later 
developed by Lazarus et al. [37]. The "edges" of the polygonal schemata produced by all 
these algorithms are (possibly overlapping) paths in the 1-skeleton of the input manifold. 
We will modify one of the algorithms of Lazarus et al. to construct short nontrivial cycles 
and cut graphs. 

Very recently, Colin de Verdiere and Lazarus consider the problem of optimizing canoni- 
cal polygonal schemata [10]. Given a canonical polygonal schema for a triangulated oriented 
manifold A4, their algorithm constructs the shortest canonical schema in the same homo- 
topy class. Surprisingly (in light of our Theorem 3.1) their algorithm runs in polynomial 
time under some mild assumptions about the input. As a byproduct, they also obtain a 
polynomial-time algorithm to construct the minimum-length simple loop homotopic to a 
given path. 

Surface parameterization is an extremely active area of research, thanks to numerous 
applications such as texture mapping, remeshing, compression, and morphing. For a sample 
of recent results, see [1, 17, 20, 48, 25, 39, 40, 45, 47, 46, 58] and references therein. In most of 
these works, surfaces of high genus are parameterized by cutting them into several (possibly 
overlapping) patches, each homeomorphic to a disk, each with a separate parameterization. 
A recent exception is the work of Gu et al. [25], which computes an initial cut graph in 
O(nlogn) time by running a shortest path algorithm on the dual of the manifold mesh, 
starting from an arbitrary seed triangle. Essentially the same algorithm was independently 
proposed by Steiner and Fischer [48]. Once a surface has been cut into a disk (or several 
disks), further (topologically trivial) cuts are usually necessary to reduce distortion [25, 45, 
47]. Many of these algorithms include heuristics to minimize the lengths of the cuts in 
addition to the distortion of the parameterization [25, 40, 45], but none with theoretical 
guarantees. 

All of our algorithms are ultimately based on Dijkstra's single-source shortest path algo- 
rithm [14, 51]. Many previous results have used Dijkstra's algorithm or one of its continuous 
generalizations [32, 41, 53] to discover interesting topological structures in 2-manifolds, such 
as cut graphs [25, 48], small handles ('topological noise') [24], texture atlases [40], contour 
trees [2, 38], and Reeb graphs [30, 48]. 
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3 Computing Minimum Cut Graphs is NP-Hard 

In this section, we prove that finding a minimum cut graph of a triangulated manifold is 
NP-hard. We consider two versions of the problem. In the weighted case, the manifold is 
assumed to be a polyhedral surface in 1R 3 and we want to compute the cut graph whose total 
Euclidean length is as small as possible. In the unweighted case, we are compute the cut 
graph with the minimum number of edges; the geometry of the manifold is ignored entirely. 

Both reductions are from the rectilinear Steiner tree problem: Given a set P of n points 
from a m x m square grid in the plane, find the shortest connected set of horizontal and 
vertical line segments that contains every point in P. This problem is NP-hard, even if 
m is bounded by a polynomial in n [23]. Our reduction uses the Hanan grid of the points, 
which is obtained by drawing horizontal and vertical lines through each point, clipped to the 
bounding box of the points. At least one rectilinear Steiner tree of the points is a subset of 
the Hanan grid [28]. 

Theorem 3.1. Computing the length of the minimum (weighted or unweighted) cut graph 
of a triangulated punctured sphere is NP-hard. 

Proof: First consider the weighted case, where the weight assigned to each edge of the 
manifold is its Euclidean length. Let P be a set of n points in the plane, with integer 
coordinates between 1 and m. We construct a punctured sphere in 0(n 2 ) time as follows. 
Assume that P lies on the xy-plane in 1R 3 . We modify the Hanan grid of P by replacing 
each terminal with a square of width l/2n, rotated 45 degrees so that its vertices lie on 
the neighboring edges. These squares will form the punctures. We then attach a basin 
under each face / of the modified Hanan grid, by joining the boundary of / to a slightly 
scaled copy of / on the plane z = — n 2 . We also attach a basin of depth n 2 + 1 to the 
boundary of the entire modified Hanan grid. The side facets of each basin are trapezoids. 
The basins are tapered so that adjacent basins intersect only on the modified Hanan grid. 
Triangulating this surface arbitrarily, we obtain a polyhedral sphere Ai with n punctures 
and overall complexity 0(n 2 ). See Figure 2. 

Let G* be a minimum weighted cut graph of Ai. We easily observe that G* contains only 
"long" edges from the modified Hanan grid and contains at least one vertex of every puncture. 
Thus, the edges of G* are in one-to-one correspondence with the edge of a rectilinear Steiner 
tree of P. 

For the unweighted case, we modify the original mx m integer grid instead of the Hanan 
grid. To create a punctured sphere, we replace each terminal point with a small diamond 
as above. We then fill in each modified grid cell with a triangulation, chosen so that the 
shortest path between any two points on the boundary of any cell stays on the boundary of 
that cell; this requires a constant number of triangles per cell. The resulting manifold Ai' has 
complexity 0{m 2 ). By induction, the shortest path between any two points on the modified 
grid lies entirely on the grid. Thus, any minimal unweighted cut graph of Ai' contains only 
edges from the modified grid. It follows that if the minimum unweighted cut graph of Ai' 
has r edges, the length of any rectilinear Steiner tree of P is exactly r. □ 
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We can easily generalize the previous proof to manifolds with higher genus, with or 
without boundary, oriented or not, by attaching small triangulated tori or cross-caps to any 
subset of punctures. 

Theorem 3.2. Computing the length of the minimum (weighted or unweighted) cut graph 
of a triangulated manifold with boundary, with any fixed genus or with any fixed number of 
boundary components, is NP-hard. 



4 Computing Minimum Cut Graphs Anyway 

We now describe an algorithm to compute the minimum cut graph of a polyhedral manifold 
in n°( g+k } time. For manifolds with constant Euler characteristic, our algorithm runs in 
polynomial time. 

Our algorithm is based on the following characterization of the minimum cut graph as 
the union of shortest paths. A branch point of a cut graph is any vertex with degree greater 
than 2. A simple path in a cut graph from one branch point or boundary point to another, 
with no branch points in its interior, is called a cut path. 

Lemma 4.1. Let M. be a polyhedral 2-manifold, possibly with boundary, and let G* be a 
minimum cut graph of M. . Any cut path in G* can be decomposed into two equal-length 
shortest paths in M. 1 . 

Proof: Let G be an arbitrary cut graph of M., and consider a cut path between two (not 
necessarily distinct) branch points a and c of G. Let b be the midpoint of this path, and let 
a and f3 denote the subpaths from b to a and from b to c, respectively. Note that b may lie 
in the interior of an edge of .Mi. Finally, suppose a is not the shortest path from b to a 
in M. i . To prove the lemma, it suffices to show that G is not the shortest cut graph of A4 . 
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Let a' be the true shortest path from b to a. Clearly, a' is not contained in G. Walking 
along a' from b to a, let s be the first vertex whose following edge is not in G, and let t be 
the first vertex in G whose preceding edge is not in G. (Note that s and t may be joined by a 
single edge in Ji4\G.) Finally, let a' C a' be the true shortest path from s to t. Equivalently, 
a' is the first maximal subpath of a' whose interior lies in M. \ G. See Figure 3. 




Figure 3. If the dashed path from a to & is shorter than a, then the cut graph can be shortened by cutting along 
a' and regluing along a. 

The subpath a' cuts M. \ G into two smaller disks. We claim that some subpath o of 
either a or ft appears on the boundary of both disks and is longer than o' . Our claim implies 
that cutting A4\G along a' and regluing the pieces along a gives us a new polygonal schema 
with smaller perimeter, and thus a new cut graph shorter than G. See Figure 3 for an 
example. 

We prove our claim by exhaustive case analysis. First consider the case where the mani- 
fold Ai is orientable. We can subdivide the entire boundary of the disk Ai\G into six paths 
labeled consecutively a, ft, 7, ft, a, 8. Here, a and ft are the corresponding copies of a and ft 
in the polygonal schema. Because M. is orientable, a and a have opposite orientations, as 
do ft and ft. Either or both of 7 and 5 could be empty. See the lower left part of Figure 3. 
The subpath a' can enter the interior of the disk M. \ G from four of these six paths (a, ft, 
a, and ft) and leave the interior of the disk through any of the six paths. 

Suppose a' enters the interior of Ai \ G from a; the other three cases are symmetric. 
Figure 4 shows the six essentially different ways for a' to leave the interior of A4 \ G. In 
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Figure 4. Six cases for the proof of Lemma 4.1 for orientable manifolds; all other cases are reflections of these. 
In each case, some subpath of a or (3 appears on the boundary of both sub-disks. 



each case, we easily verify that after cutting along a', some subpath a of either a or (3 is on 
the boundary of both disks. Specifically: 

(i) If a' leaves through a or a (see Figures 4(c) and 4(a), respectively), then a is the 
subpath of a from s to t. Since both a and a' have the same endpoints and a' is a 
shortest path, a must be longer than a'. 

(ii) If a' leaves through (3 or (3 (see Figures 4(f) and 4(d), respectively), then a is the 
subpath of (3 from b to t. Indeed, p = a[b, s] ■ a' is the shortest path from b to t, and as 
such \a'\ < \p\ < \{3[b, t]\ = \a\. (Here • denotes path concatenation, and a[x, y] denotes 
the subpath of a from x to y.) 

(iii) If a' leaves through 7 (see Figure 4(e)), then a = (3. Clearly, \a'\ < \a\ = |/3| = \a\. 

(iv) Finally, if a' leaves through 5 (see Figure 4(b)), then a is the subpath of a from a to s. 
Clearly, \a'\ < \a[s, a]\ < \a\. 

If M. is non-orientable, the path a(3 could appear either with the same orientation or with 
opposite orientations on the boundary of the disk M. \ G. If the orientations are opposite, 
the previous case analysis applies immediately. Otherwise, the boundary can be subdivided 
into six paths labeled consecutively a, (3, 7, a, (3, 5. Without loss of generality, a' enters the 
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Figure 5. Six additional cases for the proof of Lemma 4.1 for non-orientable manifolds; all other cases are 
reflections or rotations of these. In each case, some subpath of a or j3 appears on the boundary of both sub-disks. 



interior of M. \ G from a and leaves through any of these six paths. The six cases are 
illustrated in Figure 5. Again, we easily verify that in each case, some subpath a of either a 
or (3 is on the boundary of both disks. We omit further details. □ 

For any cut graph G of a manifold Ai, we define the corresponding reduced cut graph G 
as follows. First we remove any topologically trivial cuts; that is, we repeatedly remove any 
edge with a vertex of degree 1 that is not on the boundary of A4. We then augment the 
cut graph by adding all the boundary edges of AL Finally, we contract each maximal path 
through degree-2 vertices into a single edge. The resulting reduced cut graph G is 2-edge 
connected, and each of its vertices has degree at least 3. Every vertex of G is either a branch 
point or a boundary point of G, and every edge of G corresponds to either a cut path or 
a boundary path in G. However, in general, not all branch points and cut paths in G are 
represented in G. 

Lemma 4.2. Let M. be a polyhedral 2-manifold with genus g and k boundary components. 
Any reduced cut graph G of M. is connected, has between max{l, k} and 4g + 2k — 2 vertices, 
and has between g + max{0, 2k — 1} and 6g + 3 k — 3 edges. 

Proof: Let G' be the cut graph corresponding to G, after all the trivial cuts have been 
removed. The boundary of the polygonal schema M. \ G' can be partitioned into cut paths 
and boundary paths, each corresponding to an edge in G. Thus, G is connected. 
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Let v and e denote the number of vertices and edges in G, respectively. If any vertex 
in G has degree d > 4, we can replace it with d — 3 trivalent vertices and d — 3 new edges of 
length zero. Thus, in the worst case, every vertex in G has degree exactly 3, which implies 
that 3v = 2e. Since G is embedded in Ai with a single face, Euler's formula implies that 
v — e + 1 = x = 2 — 2# — A; if .M is orientable, and v — e + 1 = x = 2 — g — A; if .M is 
non-orientable. It follows that t> < Ag + 2 A; — 2 and e < 6g + 3 A; — 3, as claimed. 

On the other hand, G has at least one vertex on each boundary component of Ai, and 
at least one vertex even if At has no boundary, so v > max{l,A;}. Thus, Euler's formula 
implies that 2 — 2g — k = v — e + 1 > max{l, A:} — e + 1 if Ai is orientable, or equivalently, 
^ > 2g + max{0, 2A; — 1}. Similarly, if Ai is non-orientable, Euler's formula implies that 
e > # + max{0,2A;- 1}. □ 

Our minimum cut graph algorithm exploits Lemma 4.1 by composing potential minimum 
cut graphs out of 0(g + k) shortest paths. Unfortunately, a single pair of nodes in Ai could 
be joined by 2 n ^ shortest paths, in 2 n( - 9+k ^ different isotopy classes, in the worst case. To 
avoid this combinatorial explosion, we can add a random infinitesimal weight e-w(e) to each 
edge e. The Isolation Lemma of Mulmuley, Vazirani, and Vazirani [42] implies that if the 
weights w(e) are chosen independently and uniformly from the integer set {1, 2, . . . , n 2 }, all 
shortest paths are unique with probability at least 1 — 1/n; see also [9, 34]. 4 

We are now finally ready to describe our minimum cut graph algorithm. 

Theorem 4.3. The minimum cut graph of a polyhedral 2-manifold Ai with genus g and k 
boundary components can be computed in time n 0( - 9+k K 

Proof: We begin by computing the shortest path between every pair of vertices in Ai 
in 0(n 2 log n) time by running Dijkstra's single-source shortest path algorithm for each 
vertex [14, 31], breaking ties using random infinitesimal weights as described above. Once 
these shortest paths and midpoints have been computed, our algorithm enumerates by brute 
force every possible cut graph that satisfies Lemmas 4.1 and 4.2, and returns the smallest 
such graph. 

Each cut graph is specified by a set V of up to Ag + 2 A: — 2 vertices of Ai, a set E of up to 
6(7 + 3A: — 3 edges of Ai, a trivalent multigraph G with vertices V , and a assignment of edges 
in E to edges in G. Each edge (v, w) of G is assigned a unique edge e G E to define the 
corresponding cut path in Ai. This cut path is the concatenation of the shortest path from 
v to e, e itself, and the shortest path from e to w. If the midpoint of this cut path is not 
in the interior of e, we declare the cut path invalid, since it violates Lemma 4.1. (Because 
shortest paths between vertices are unique, the midpoint of any cut path in the minimal cut 
graph must lie in the interior of an edge.) If all the cut paths are valid, we then check that 

4 Alternately, if we choose w(e) uniformly from the real interval [0,1], shortest paths are unique with 
probability 1. This may sound unreasonable, but recall that no polynomial-time algorithm is known to 
compare sums of square roots of integers in any model of computation that does not include square root 
as a primitive operation [7]. Thus, to compute Euclidean shortest paths in a geometric graph with integer 
vertex coordinates, we must either assume exact real arithmetic or (grudgingly) accept some approximation 
error [22]. 
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Figure 6. From left to right: non-separating, essential but separating, and trivial cycles on a 2-manifold. 

every pair of cut paths is disjoint, except possibly at their endpoints, and that removing all 
the cut paths from M. leaves a topological disk. 

Our brute-force algorithm considers 0(n i9+2k ~ 2 ) different vertex sets V, 0(n 6g+3k ~ 2 ) 



Thus, n°^ 9+k ^ potential cut graphs are considered altogether. The validity of each potential 



5 Finding Short Nontrivial Cycles 

As a step towards efficiently computing approximate minimum cut graphs, we develop al- 
gorithms to compute shortest and nearly-shortest nontrivial cycles in arbitrary 2-manifolds, 
possibly with boundary. Although our most efficient approximation algorithm for cut graphs 
requires only approximate shortest nontrivial cycles on manifolds without boundary, we be- 
lieve these algorithms are of independent interest. 

We distinguish between two type of nontrivial simple cycles. A simple cycle 7 in M. 
is non-separating if M. \ 7 has only one connected component. A simple cycle 7 in M. is 
essential if it is not contractible to a point or a single boundary cycle of M.. Every non- 
separating cycle is essential, but the converse is not true. Formally, non-separating cycles 
are homologically nontrivial, and essential cycles are homotopically nontrivial. See Figure 6. 

As can be seen immediately from Figure 6, it is not possible to determine whether a cycle 
is non-separating, essential, or trivial by examining only a local neighborhood. Dey and 
Schipper [13] describe an algorithm to determine whether an arbitrary circuit is contractible 
in 0(n) time; their algorithm begins by computing an arbitrary cut graph of the manifold. 
Fortunately, as we shall see shortly, we can simplify their algorithm considerably when the 
given cycle is simple. 

5.1 Shortest Cycles 

We begin by describing how to find the shortest nontrivial cycle through a given vertex. 
Our algorithm uses a combination of Dijkstra's single-source shortest path algorithm [14] 
and a modification of the canonical polygonal schema algorithm of Lazarus et al. [37]. Our 





□ 
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algorithm is similar to the approach taken by Guskov and Wood [24] to find and remove 
small handles from geometric models reconstructed from noisy data; see also [57]. 

The algorithm of Lazarus et al. builds a connected subset S of a triangulated manifold 
without boundary, starting with a single triangle and adding new triangles on the boundary of 
S one at a time. If a new triangle intersects the boundary of S in more than one component, 
the algorithm checks which of the following three cases holds: (1) At \ S is connected; 
(2) neither component of At \ S is a disk; or (3) one component of At \ S is a disk. In the 
final case, the algorithm adds the disk component to S and continues searching the other 
component of At \ S. If we run this algorithm until either case (1) or case (2) holds, the 
total running time is 0(n). See Lazarus et al. [37] for further details. 

First we describe a straightforward generalization of the algorithm used by Lazarus to 
determine the structure of At \ S. 

Lemma 5.1. Let At be a connected polyhedral 2-manifold At whose genus g and number 
of boundary components k are known. Any simple cycle 7 in At\ can be classified as non- 
separating, essential but separating, or trivial in 0(m) time, where m is the complexity of 
the smaller component of At \ 7. fn particular, if 7 is non-separating, the running time is 
0(n). 

Proof: We perform two simultaneous depth-first searches, starting from either side of any 
vertex of 7. If either search meets a vertex already visited by the other, 7 is a non-separating 
cycle. The running time in this case is trivially 0(n). 

Conversely, if one search halts without meeting the other, 7 is a separating cycle. Let C 
be the smaller component of At \ 7. To determine whether 7 is essential, we compute the 
Euler characteristic and the number of boundary components of C, which we denote by x(C) 
and k(C), respectively. (This can be done on the fly during the depth-first search phase.) 
Let x denote the Euler characteristic of At. The cycle 7 is contractible if and only if one of 
the following conditions holds. 

• C is a disk, or equivalently, x(C) = 1. 

• C is an annulus, or equivalently, x(C) = and k(C) = 2. 

• At \ C is a disk, or equivalently, x(C) = x ~ 1- The equivalence follows from the 
inclusion-exclusion formula x{A U B) — x(A) + x{B) — x{A n B). 

• M. \ C is an annulus, or equivalently, — x(-M) and k(C) = k(M). 

If none of these conditions hold, then 7 is essential. The total running time if 7 is a separating 
cycle is 0(m), where m is the complexity of C. □ 

Note that we can simplify this algorithm slightly if the manifold has no boundary, since 
in that case neither component of Ai \ 7 is an annulus. We will use this simplification in our 
approximate cut graph algorithm. 

Lemma 5.2. Let u be a vertex of a polyhedral 2-manifold Ai, possibly with boundary. The 
shortest essential cycle in Ai\ that contains u can be computed in O(nlogn) time. 
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Proof: We find the shortest essential cycle through u by simulating a circular wave expand- 
ing from u. Whenever the wave touches itself, either we have the shortest essential cycle 
through u, or one component of the wave bounds a disk in Ai and we can continue expanding 
the other component. 

We modify the algorithm of Lazarus et al. [37] as follows. First, S is no longer a set of 
triangles but a more general connected subset of vertices, edges, and facets of M. Initially, 
S contains only the source vertex u. Second, we use Dijkstra's algorithm to determine the 
order for edges to be added. We add a facet to S only when all its vertices have been added 
to S, either directly or as part of another facet. We run the algorithm from Lemma 5.1 
whenever S is no longer simply connected, that is, when we add a new edge vw with both 
endpoints on the boundary of S. If M.\S is disconnected, we continue only if one component 
of M, \ S is a disk or an annulus, as can be checked using Lemma 5.1. In that case, we 
add the disk or annulus component of A4 \ S to S, discard the vertices of that component 
from the Dijkstra priority queue, and continue searching in the other component. If M. \ S 
is connected, or if neither component is a disk or an annulus, we have found the shortest 
essential cycle through u, consisting of the shortest path from u to v, the edge vw, and the 
shortest path from w to u. 

Dijkstra's algorithm requires 0(n log n) time. Each time we find a trivial cycle, we spend 
0(m) time and discard a disk with complexity at least m. Thus, the total time spent 
performing cycle classification and maintaining the wavefront set S is 0{n). Thus, the total 
running time of our algorithm is O(nlogn). □ 

Running this algorithm once for every vertex of M. immediately gives us the shortest 
essential cycle. 

Corollary 5.3. Let M. be a polyhedral 2-manifold, possibly with boundary. The shortest 
essential cycle in M.\ can be computed in 0(n 2 logn) time. 

A simple modification of our algorithm allows us to find shortest non-separating cycles 
in the same asymptotic time. 

Lemma 5.4. Let u be a vertex of a polyhedral 2-manifold A4, possibly with boundary. The 
shortest non-separating cycle in M. 1 that contains u can be computed in 0(n\ogn) time. 

Proof: The only change from the previous algorithm is that if we discover an essential 
separating cycle, we continue recursively in both components of M\S. The cost of Dijkstra's 
algorithm is still O(nlogn), but we now must spend extra time in the cycle-classification 
algorithm of Lemma 5.1. As before, the total time spent finding trivial cycles is 0(n), since 
we can charge the search time to the discarded components. 

Let T(n,g) denote the total time spent finding separating essential cycles. This function 
satisfies the recurrence 



T(n, g) < T(m, h) + T(n - m, g - h) + 0(m), 
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where m < n/2 is the complexity of the smaller component of M. \ S and h is its genus. The 
base case of the recurrence is T(n, 1) = 0, since every essential cycle on a genus-1 surface is 
non-separating. 

Similar recurrences appear in the analysis of output-sensitive planar convex hull algo- 
rithms [5, 8, 33, 56], suggesting that the solution to our recurrence is T(n,g) = 0(n\ogg). 
Indeed, we can prove this by induction as follows. Suppose 

T(n, g) < T(m, h) + T{n — m, g — h) + cm 

for some constant c. We claim that T(n,g) < cnlgg. The inductive hypothesis implies that 

T(n, g) < cm lg h + c(n — m) \g{g — h) + cm 

< max (cm lg h + c(n — m) lg{g — h)) + cm. 

A simple application of derivatives implies that the right hand side of this inequality is 
maximized when h = mg/n. Thus, 

ml X 1 m 9 I \ 1 ( n ~ m )9 

Tin, g) < cm lg h cm — m) lg h cm 

n n 

= cnlgg + cm lg m + c(n — m) lg(n — m) — cn lg n + cm. 
Since m < n/2 and n — m < n, we can simplify this inequality to 

T{n,g) < cn In g + cm lg(n/2) + c(n — m) \gn — cnXgn + cm = cnlgg, 
completing the proof. 

Thus, the total time spent in the cycle-classification phase of our algorithm is 0(n\ogg). 
Since g < n, this is dominated by the cost of maintaining the Dijkstra priority queue. □ 

Corollary 5.5. Let M. be an polyhedral 2-manifold, possibly with boundary. The shortest 
non-separating cycle in A4i can be computed in 0(n 2 \ogn) time. 

5.2 Nearly-Shortest Cycles 

As we will argue in the next section, computing short nontrivial cycles is the bottleneck in 
our approximate cut graph algorithm. Fortunately, exact minimum cycles are not necessary 
for our results. We can speed up our cut graph algorithm, without significantly increasing 
the approximation factor, by searching for a nontrivial cycle at most twice as long as the 
shortest. Our approximation algorithm assumes that the manifold M. has no boundary; 
fortunately, as we shall see in the next section, this is sufficient for our purposes. 

Our approximation algorithm works as follows. First, we compute a set of shortest paths 
(in fact, a cut graph) that intersects every non-separating cycle in the manifold M.. Then 
we contract each shortest path n in this set to a point and find the shortest nontrivial cycle 
through that point, as described by Lemmas 5.2 and 5.4. 
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Lemma 5.6. Let ir be a shortest path between two vertices in a polyhedral 2-manifold Ai, 
and let 7* be the shortest essential (resp. non- separating) cycle in Aii that intersects tt. In 
0(n\ogn) time, one can compute an essential (resp. non-separating) cycle 7 in Ai such that 
|7| < 2|7*|. 

Proof: Let At' be the manifold obtained by contracting the shortest path ix to a single 
vertex v. Because tt has no cycles, Ai' has the same topological type as Ai. Let 7' be 
the shortest essential (resp. non-separating) cycle in Ai' that passes through v. Clearly, 
It'I < |t*|- We can compute this cycle in O(nlogn) time by Lemma 5.2 (resp. Lemma 5.4). 

We construct a cycle 7 in Ai by concatenating two paths a and f3, where a contains the 
edges of 7' and (3 is the subpath of n between the endpoints of a. The sequence of edge 
contractions that transforms Ai to Ai' also transforms 7' to 7. Hence, 7' is an essential cycle 
of Ai. Because j3 is a subpath of a shortest path, f3 is actually the shortest path between 
the endpoints of a, so \(3\ < \a\ = \j'\. It follows that (7] = |a| + \(3\ < 2\j'\ < 2\j*\. □ 

This lemma suggests a natural algorithm for finding a short nontrivial cycle: Compute 
a set of shortest paths that intersect every non-separating cycle (and thus every essential 
cycle), and then run the algorithm from Lemma 5.6 for each path in this set. 

Lemma 5.7. Let M be a polyhedral 2-manifold without boundary. In 0(n\ogn) time, one 
can compute a set IT of 0(g) shortest paths on Aii such that every non-separating cycle 
(and thus every essential cycle) in Aii intersects at least one path in LT. 

Proof: We compute a cut graph G as follows. First we compute a shortest-path tree T 
from an arbitrary initial vertex v using Dijkstra's algorithm. We then compute an arbitrary 
spanning tree T* of the dual of Ai \ T, that is, the graph whose vertices are facets of At 
and whose edges join pairs of facets that share a common edge not in T. Analysis similar to 
Lemma 4.2 implies that there are 0(g) edges that do not appear in T and whose dual edges 
to not appear in T*. Call this set of unclaimed edges E. Let LT be the set of 0(g) shortest 
paths from v to the endpoints of E; these paths are all in T. Finally, let G = LT U E'. 

We easily observe that Ai \ G is a topological disk, so G is a cut graph. It follows that 
every non-separating cycle in Ai\ intersects G. Since every vertex of G is also a vertex of 
some path in LT, every non-separating cycle in Aii intersects at least one path in LT. □ 

Notice that this algorithm does not work if Ai has a boundary, since the dual graph of 
Ai \ T could be disconnected. 

Corollary 5.8. Let Ai be a polyhedral 2-manifold with genus g and no boundary, and let 7* 
be its shortest essential (resp. non-separating) cycle. In O(gnlogn) time, one can compute 
an essential (resp. non-separating) cycle 7 in Ai\ such that \y\ < 2|7*|. 

Proof: We construct a set LT of 0(g) shortest paths, at least one of which is guaranteed to 
intersect 7*, as described in the previous lemma. Then for each path tt e II, we contract 7r 
to a point and find the shortest nontrivial cycle through that point in O(nlogn) time, as 
described by Lemma 5.6. □ 
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6 Approximate Minimum Cut Graphs 

We now describe a simple polynomial-time greedy algorithm to construct an approximate 
minimum cut graph for any polyhedral manifold Ai. 

To handle manifolds with boundary, it will be convenient to consider the following simpli- 
fied form. Given a manifold Ai with genus g and k boundary components, the corresponding 
punctured manifold (Ai, P) consists of a manifold Ai with the same genus as Ai but without 
boundary, and a set P of k points in Ai, called punctures. To construct Ai, we contract 
every boundary component of Ai to a single point, which becomes one of the punctures 
in P. 5 If any vertex of Ai has multiple edges to the same boundary component, Ai contains 
only the edge with smallest weight, breaking ties using the Isolation Lemma as above. If Ai 
has no boundary, then Ai — Ai and P = 0. 

Our goal now is to compute the minimum cut graph of Ai that touches every puncture 
in P; henceforth, we call this simply the minimum cut graph of (Ai,P). This reduction is 
motivated by the following trivial observation. 

Lemma 6.1. The minimum cut graph of any polyhedral 2-manifold Ai has the same length 
as the minimum cut graph of (Ai, P). 

Our approximation algorithm works as follows. We repeatedly cut along short nontrivial 
cycles until our surface becomes a collection of punctured spheres, connect the punctures on 
each component by cutting along a minimum spanning tree, and finally (if necessary) reglue 
some previously cut edges to obtain a single disk. The resulting cut graph is composed of a 
subset of the edges of the short nontrivial cycles and all the edges of the minimum spanning 
forest. 

6.1 Using Short Nontrivial Cycles 

The first component of our algorithm is a subroutine to compute approximately shortest 
nontrivial cycles, described by Corollary 5.8. As required by that algorithm, the input 
manifold Ai has no boundary; the punctures are completely ignored. 

The following argument relates the length of the shortest nontrivial cycles to the length 
of the minimum cut graph. 

Lemma 6.2. Let G be any cut graph of a polyhedral 2-manifold Ai with genus g and no 
boundary. The shortest cycle in G contains 0((log g)/g) of the total length of G. 

Proof: First consider the reduced cut graph G, constructed by repeatedly contracting any 
edge with a vertex of degree less than three, as in Section 4. Every vertex in G has degree 
at least 3. Without loss of generality, assume that every vertex in G has degree exactly 3, 
splitting each high-degree vertex into a tree of degree-3 vertices if necessary, as in the proof 
of Lemma 4.2. A straightforward counting argument implies that any trivalent graph whose 

5 We could simulate this contraction by artificially assigning every boundary edge of AI a weight of zero, 
although this would require a few simple changes in our algorithms. 
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girth (minimum cycle length) is c must have at least 3\/2 • 2 C / 2 — 2 vertices if c is odd, and 
at least 2 • 2 C//2 — 2 vertices if c is even [6]. By Lemma 4.2, G has at most 4g — 2 vertices, so 
G must have a cycle 7 with most 2(lgg + 1) = O(logg) edges. 

Starting with G = G, we inductively define a sequence of reduced graphs G\, G 2 , ■ ■ ■ as 
follows. For each % > 0, let % denote the shortest cycle in We obtain Gi by reducing 

the graph Gj_i \ or equivalently, removing the vertices of % and all their edges, and 
then contracting \^\ nearby length-2 paths to single edges. Our earlier argument implies 
that each cycle 7$ has at most 2(lgg + 1) edges. Thus, for each i, we have \E{Gj)\ = 
\E(Gi-i)\ — 6(lgg + 1). Lemma 4.2 implies that the original reduced cut graph G has at 
least g edges, so we can repeat this process at least g/6(\gg + 1) times. 

Let 7j denote the cycle in the original cut graph G corresponding to G^. By our con- 
struction, 7j and 7, are disjoint for all i 7^ j, so we have a set of at least g/6(lgg + l) disjoint 
cycles in G. At least one of these cycles has length at most 6(lgg + l)/g = 0((log g)/g) 
times the total length of G. □ 

Since every cycle in the minimum cut graph is non-separating, and therefore essential, 
we immediately have the following corollary. 

Corollary 6.3. For any polyhedral 2-manifold M. with genus g and no boundary, both the 
length of the shortest non-separating cycle and the length of the shortest essential cycle are 
at most 0((log g)/g) times the length of the minimum cut graph of M.. 

6.2 Puncture- Spanning Trees 

The second component of our cut graph algorithm is a subroutine to compute the minimum 
puncture- spanning tree of a punctured manifold (Ai, P), that is, the minimum spanning tree 
of the punctures P in the shortest-path metric of M.\. 

Lemma 6.4. The minimum puncture-spanning tree of any punctured polyhedral 2-manifold 
(Ai,P) can be computed in 0(n\ogn) time. 

Proof: We simulate Kruskal's minimum spanning tree algorithm [36, 51] by adding shortest 
puncture-to-puncture paths one at a time, in increasing order of length. To compute the 
shortest paths, we simultaneously propagate wavefronts from all k punctures using Dijkstra's 
algorithm. Whenever two wavefronts (i.e., two growing shortest-path trees) collide, we 
add a new edge to the evolving minimum spanning tree and merge those two wavefronts. 
To implement this algorithm efficiently, we maintain the wavefronts in a union-find data 
structure. The resulting running time is 0{n log n). □ 

This is essentially the algorithm proposed by Takahashi and Matsuyama [50] to compute 
approximate Steiner trees in arbitrary graphs. The same algorithm was also recently used 
by Sheffer [45] and by Levy et al. [40] to compute cut graphs, where surface features with 
high discrete curvature play the role of punctures. 
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Lemma 6.5. The length of the minimum puncture-spanning tree of any punctured polyhe- 
dral 2-manifold (M, P) is at most twice the length of any cut graph of (M, P). 

Proof: The minimum Steiner tree of P is the subgraph of Aii of minimum total weight 
that includes every point in P. Since any cut graph of (Ai,P) must touch every puncture, 
no cut graph is shorter than this minimum Steiner tree. On the other hand, the minimum 
spanning tree of P has at most twice the length of the minimum Steiner tree [35, 50]. □ 

6.3 Analysis 

We now have all the components of our greedy cut graph algorithm. At any stage of the 
algorithm, we have a punctured manifold (M.,P). Our algorithm repeatedly cuts along 
a short non-separating cycle of Ai, using Corollary 5.8. This cut creates one or two new 
boundary circles, which we collapse to new punctures. When the manifold is reduced to a 
collection of punctured spheres, we cut along the minimum puncture-spanning tree of each 
component using the algorithm in Lemma 6.4. 

Each non-separating cycle cut reduces the genus of Ai by 1. This immediately implies 
that our algorithm performs exactly g cycle cuts, so the overall running time is 

g ■ 0(gn\ogn) + O(nlogn) = 0(g 2 n\ogn). 

For any graph X, let \X\ denote its total length. Let G* denote the minimum cut graph 
of (A4,P). Let (Mi, Pi) denote the punctured manifold after g — i cycle cuts have been 
performed, so Aii has genus i, and let G* denote the minimum cut graph of Aii, ignoring 
the punctures Pj. Since collapsing edges cannot increase the minimum cut graph length, we 
have < \G*\ < \G*\ for all i. 

Let 7, denote the short non-separating cycle of Aii found by Corollary 5.8. (We easily 
observe that any cut graph of Aii must intersect this cycle.) Lemma 6.2 and Corollary 5.8 
imply that |7j| < 0((\ogi)/i) ■ \G*\ for all i and j. Summing over all g cuts, we conclude 
that the total length of all cycle cuts is at most 

J20((\ogt)/t)-\G*\=0(\og 2 g)-\G*\. 

i=i 

Similarly, Lemma 6.5 implies that the minimum puncture-spanning forest has length at 
most 2|G*|. Finally, regluing previously cut edges to obtain a single disk only reduces the 
length of the final cut graph. Thus, the final cut graph computed by our algorithm has 
length at most 0(log 2 g) ■ \G*\. 

Theorem 6.6. Given a polyhedral 2-manifold At with genus g and k boundary components, 
an O (log 2 g) -approximation of its minimum cut graph can he constructed in 0(g 2 n\ogn) 
time. 

Cutting along short essential cycles instead of short non-separating cycles leads to exactly 
the same asymptotic running time and approximation bounds, although the algorithm and 
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its analysis are slightly more complicated. For purposes of analysis, we can divide the 
algorithm into phases, where in the ith phase, we cut along a short essential cycle of every 
component of the manifold that has genus i. Essential cycle cuts can separate the manifold 
into multiple components, but since each component must have nontrivial topology, the 
algorithm performs at most g — 1 separating cuts. At the end of the algorithm, if necessary, 
we reglue along some previously cut edges to obtain a single topological disk. We refer to 
the earlier version of this paper for further details [19]. 



7 Open Problems 

We have developed new algorithms to compute exact and approximate minimal cut graphs 
for manifold surfaces with arbitrary genus and arbitrary boundary complexity. Our approx- 
imation algorithm is particularly simple. 

Our results suggest several open problems, the most obvious of which is to improve the 
running times and approximation factors of our algorithms. Is the minimum cut graph 
problem fixed-parameter tractable [15]? That is, can we compute exact minimum cut graphs 
in time f(g,k)- for some function /? The similarity to the Steiner problem offers some 
hope here, since the minimum Steiner tree of k nodes in an n-node graph can be computed 
in 0(3 k n + 2 k n 2 + n 3 ) time [16, 27]. 

The approximation algorithm of Theorem 6.6 is somewhat indirect. It computes a short 
cut graph by repeatedly computing a 'reasonable' cut graph and then extracting a short 
nontrivial cycle that interacts with this cut graph. It is natural to conjecture that one can 
compute such a short cut graph directly, resulting in a faster algorithm. In particular, we 
conjecture that an approximately minimum cut graph can be computed in O(gnlogn) time. 

How well can we approximate the minimum cut graph in nearly-linear time? There are 
several simple heuristics to compute 'good' cut graphs in O(nlogn) time, such as the dual 
shortest-path algorithm used by Gu et al. [25] and by Steiner and Fischer [48], and the 
algorithm described in the proof of Lemma 5.7. How well do these algorithms approximate 
the minimum cut graph? 

More generally, is there a simple, practical, 0(l)-approximation algorithm, like the min- 
imum spanning tree approximation of Steiner trees? In fact, it might be that our algorithm 
provides such an approximation, as our current analysis seems to be far from tight. Un- 
fortunately, the general Steiner tree problem is MAXSNP-hard [4], so an efficient (1 + e)- 
approximation algorithm for arbitrary e > seems unlikely. 

Several authors have pointed out apparent tradeoffs between the quality of parameteri- 
zations and the length of the required cut graph; see, for example, Sorkine et al. [47]. How 
hard is it to compute the (approximately) shortest cut graph required for a parameterization 
whose distortion is less than some given limit? Conversely, how hard is it to (approximately) 
minimize the distortion of a parameterization, given an upper bound on the permitted length 
of the surface cuts? The complexity of these problems clearly depends the which distortion 
measure is used, but we expect almost any variant of this problem to be NP-hard. 
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Finally, can our ideas be applied to other useful families of curves on manifolds, such as 
homology generators (families of 2g cycles that intersect in g pairs) and pants decompositions 
(maximal sets of pairwise disjoint essential cycles [52])? 
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